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ABSTRACT 

■Sf&l 

In many recent papers concerned with providing an expla- 
nation for the geomagnetic anomaly, agreement with measured 
data has been obtained from the equations of motion for elec- 
trons and ions when used with an empirical boundary condition, 
whereas poor agreement has resulted from attempts to numeri- 
cally integrate the commonly employed form of the continuity 
equation. We have been able to explain this discrepancy by 
demonstrating that the assumptions used to derive this for™ of 
the continuity equations do not agree with observation. 

Since the equations of motion do provide a favorable 
description for the geomagnetic anomaly, we have studied the 
possible physical models leading to the form of the equations 
used, and found that although field aligned diffusive equilib- 
rium provides the correct form, a more reasonable assumption 
concerning electron and ion collisions with neutrals also leads 
to the same result. We have then been able to provide a more 
realistic theoretical description of the geomagnetic anomaly 
by employing an analytic form for the boundary condition which 
is in better agreement with measurement than those previously 
used* 

Finally, by combining the equations of motion for neutrals 
electrons and ions, we have been able to indicate geomagnetic 
control for the neutral atmosphere in the lower F region of the 
ionosphere, although the exact shape of this distribution is 
unknown. 


INTRODUCTION 


In recent months, it has become increasingly evident that 
some confusion exists in the understanding of the basic physical 
mechanisms governing diffusion and the existence of the geomag- 
netic anomaly in the ionosphere. This apparent confusion arises 
by comparison of the work of Chandra (1964) , (to be referred to 
as C-I), in which it is shown that the assumption of ambipolar 
diffusion along a field line cannot lead to geomagnetic control 
of the charged particles, and such papers as Goldberg and 
Schmerling (1962, 1963), (to be referred to as GS-I and GS-II), 
and Goldberg, Kendall, and Schmerling (1964), (to be referred 
to as GKS) , in which this process does appear to produce geo- 
magnetic control of the charged particle density in the iono- 
sphere . 

The purpose of this paper is to describe and resolve the 
confusion which exists in the field at this moment, and then to 
point out the new problems with which we must contend in order 
to derive and apply the diffusion equation to ionospheric prob- 
lems correctly. In addition, a section will be devoted to an 
improved theoretical description of the geomagnetic anomaly by 
using an analytic expression for the vertical electron density 
distribution at the equator which is more in accordance with 
measurement than the simple Chapman type distribution employed 
in GKS. 

FUNDAMENTAL EQUATIONS AND DEFINITIONS 

The major cause of confusion appears to lie in the appli- 
cation of two phrases, viz. ambipolar diffusion and diffusive 
equilibrium . Let us investigate and discuss each of these terms 
to determine how loose usage of them has led to some of the 
current problems of misunderstanding. 

In the normal sense, ambipolar diffusion refers to a plasma 
in which the negative (electrons) and positive (ions) charges do 
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not move independently due to the influence of the electric 
field caused by their Coulomb interactions. In this medium, 
the electrons and ions drift in pairs and this motion of 
electron-ion pairs is referred to as ambipolar diffusion. The 
condition for ambipolar diffusion in a neutral plasma is thereby 

V e = v A = V (1) 

where v is macroscopic velocity and the subscripts e and i 
refer to electrons and ions respectively. When 

v = 0 (2) 

the condition for diffusive equilibrium is satisfied. 

The implications of (1) are quite straightforward, as shown 
in C-I. In an isothermal atmosphere and in the presence of a 
magnetic field, this requires Vx(vx§) =0. In particular, the 
assumption of field aligned plasma diffusion (vx§ = 0) can only 
be satisfied for the trivial case, v = 0, resulting in a hydro- 
static distribution of electron density independent of geomag- 
netic latitude. 

On the other hand, favorable comparison between Alouette 
topside sounder measurements and theoretical calculations of 
the geomagnetic anomaly has been obtained in GKS by assuming 
conditions of ambipolar diffusion and diffusive equilibrium 
along field lines, thereby indicating a possible conflict with 
the results in C-I. The problem resolves itself once one investi- 
gates the meaning of ambipolar diffusion and diffusive equilib- 
rium in the GKS sense. 

Let us first write the general equations of motion for 
neutrals, electrons and ions, respectively, where the subscript 
n refers to neutrals. Following C-I: 
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n.m.m 
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e l 


v ei ( W + / l ° v ln <VV " - 7p i + n i m i 5 

l n 


+e n i (2 + ^ x 8) (5) 


where n is number density, v. . is the collision frequenfey between 
the k and l x particle, m is mass, p is pressure, g is grav- 
itational acceleration, e is the absolute value of electron 
charge, 2 is electric field, and § is magnetic field. In writing 
equations (3) - (5) it is assumed that ■ — = — — . 

In the following we assume that the plasma is in a quasi- 
neutral state 


n e w n i = N (6) 

and that the electrons, ions and neutrals obey the ideal gas 
law in the ionosphere, 


p . = n .k T. 
J J J 


( 7 ) 


4 


where k is Boltzmann's constant and T is temperature. Further- 
more, we assume thermal equilibrium, i.e. 

T e = T i = T (8) 

Then, 

p e = p i = p (9) 

In addition, we assume for simplicity that 


v n « 0 (10) 

Then summation of (4) and (5) provides 
m m N m m.N 

Ven ? + _2_i v in v. = -2Vp + N(m e +m.)g + 3 x S (11) 

m e + m n + m n 

where 3 is current density, defined as 


3 = Ne (v. - v • ) 

X v 


( 12 ) 


Since we are investigating ambipolar diffusion and diffusive 
equilibrium in the GKS sense, it is desirable to write this 
equation in component form along a field line as 


m m 
n e 

in + m„ 
e n 


v . h + 
en e 


m m. 
n i 


m. + m 
i n 


; in ^ . H = [ 2k y - 7 ~ + (m e +m jL )g S (13) 


where h is a unit vector in the direction of the magnetic field. 
Let us write (13) in more familiar form by using 


- 5 - 

m e <<: m i » m n (14) 

and defining the scale height of the ionizable constituent as 
EL , where 


„ _ kT 

H i m i g 


Also, for convenience, we make the approximation 


m Ri m . 

n x 


(15) 


(16) 


Then 


-V s \ • s + r v in 7 i • B “ - kT< 


7N , r x t. 

W + 2H7' ,h 


(17) 


Finally, we write 


B = -(x sinl + x. cos I) (18) 

I* u 

where x and x„ are unit vectors in the r and 0 directions and 
r 9 

I is the magnetic dip angle, reckoned positive when the north 
seeking pole of the needle points downward. Now, if we treat 
ambipolar diffusion in the GKS sense, we simply imply that the 
electron and ion velocity components in the field direction are 
equa 1 , i . e . 


v . h = v. . h = v 


(19) 


1 1 


Applying (18) and (19) in (17), we obtain 


v = 
i i 


-kT 


.1 3N 


— i <#?# + -nr 


) + 


cos I BN 

~ nF~ ae 


] 


( 20 ) 
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where 


MV 


m v m. v . 
e en j. 1 in 

T 1 A 1 1 


Assuming that 


m v « m.v. 

e en i in 


because of (14) , we may write 


uv 


m i v in 


( 21 ) 


( 22 ) 


(23) 


Equation (20) is a familiar result derived in such papers as 
Kendall (1962) and GS-II. However, it is clearly not the result 
of ambipolar diffusion, which is given by (1), but instead, the 
result of a statement concerning the field line components of 
electron and ion velocities given by (19) . 

If we now demand 


v =0 (24) 

i i 

which is the statement implying diffusive equilibrium along a 
field line in the GKS sense, we obtain the familiar equation 

sin 1 ( N hr + mr* + Nr a? " 0 (25) 

which can also be written as 


1 dN 
N dr 



( 26 ) 
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provided we recognize that r and 0 are not independent in (25) 
but related by the dipole field conditions 

2 

r = r Q sin 0 (27) 


and 


tan 1=2 cot 0 


(28) 


It is evident that (26) can only be treated in total deriv- 
ative form if the integration is carried out along the field 
line . 

Statements concerning the components of vectors in a par- 
ticular direction, such as (19) , do not imply any conditions 
on the total vector. As a result, (25) has not required the 
assumption of any restrictions on the behavior of the velocity 
components normal to the field lines. 

Equation (26) has been the basis of describing geomagnetic 
control in the upper F-region in GKS paper. Although (26) 
has been derived assuming diffusive equilibrium along a field 
line, it is undersirable to apply this concept because it is of 
purely hypothetical nature. We now investigate other assumptions 
to find a more realistic justification for (26). 

Let us rewrite (17) as 


m v m. v . 

-TT 5 v e + ^ v i = - M < 


VN , J 

n~ + 2ir 


).H 


(29) 


i i 


1 1 


We find two ways for the right hand side of (29) toapproach zero 
The first imposes a new condition on the velocities, viz: 
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or 


m. v . v. 
V — — 1 

e 2m e v en 


V 

e 

1 i 


m. v . v . 
x in i 


-XX- 


2m v 
e en 


(30) 


a result which, although possible, would require a very special 
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condition that the electron velocity be of the order of 10 
times greater than and in the opposite direction of the ion 
velocity. 

However, if we can demand that the collision frequencies 
between electrons and neutrals and between ions and neutrals 
be sufficiently small so that the drag forces arising due to 
collision be negligible as compared to the pressure gradient, 
gravity and Lorentz forces, it is possible to derive equation 
(26) without imposing any restriction on the velocities of elec- 
trons and ions. We believe that this assumption is more realistic 
in the upper F-region where the gyro-frequencies of electrons 
and ions are much greater than their corresponding collision 
frequencies. 

Although the collision frequency assumption is physically 

more desirable, it prevents us from obtaining a simple expression 

for v e or V-; . Instead, we must return to the original 

nil _ _ 

equations of motion, (4) and (5), and solve for v and v. 

v JL 

explicitly, as has been carried out in the appendix in C-I. 

Unfortunately this introduces a very serious complication in 

the work because of the difficulty in eliminating electric field 

from the expressions of v and v. without making specific as- 

^ ^ 

sumptions about the relationship between v g and v i> The impli- 
cations of these assumptions will be discussed in the latter 
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part of this paper. In the following section we proceed to 
discuss the physical implications of equation (26) . 

THE ELECTRON DENSITY DISTRIBUTION WITH A VARIABLE 
SCALE HEIGHT 

Equation (26) can be integrated along a field line to provide 
the general solution 

o 

r cot 0 

2TT — 

N(r ,0) = f(r 0 ,TT/2)e 1 (31) 

However, if we treat T and m. constant but recognize that g 

2 1 2 
is proportional to 1/r , H.^ is then proportional to r , and we 

obtain 

r cos^ 0 
2rf.Tr) 

N(r , 0) = f(r o ,TT/2)e 1 (32) 

In both cases, f(r Q ,n/2) is an arbitrary function of height 
at the equator which cannot be determined by the equations of 
motion from which (31) or (32) are derived. The function 
f(r o>TT /2) must therefore be given as a boundary condition in 
this problem and can only be determined empirically or by use 
of additional equations governing the physics of the problem. 

Since (31) or (32) depend exclusively upon the equations 
of motion, it appears that an additional equation, sqch as the 
continuity equation, should lead to the desired boundary con- 
dition. Unfortunately, as we will show in the next section, 
the derivation and solution of the continuity equation depend 
upon a knowledge of 2 . Thus, the complexity of the problem 
becomes quite formidable and it is difficult to anticipate a 
simple method of solution at this time. 
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Instead we depend upon an empirical type boundary condi- 
tion, which may very well be the solution of the correct conti- 
nuity equation, to derive the explicit form of the electron 
density distribution. 

The incorporation of a Chapman distribution for the boundary 
condition in (31) leads to the results obtained in GKS. Since 
such a boundary condition can only be considered as a rough 
approximation to the shape of the actual vertical electron density 
distribution at the equator, it is desirable to employ an analytic 
boundary condition which more closely resembles the true height 
profiles. Chandra (1963) has proposed a modified form of the 
Chapman function which includes the effect of variable scale 
height and which is found to fit the measured vertical distri- 
bution for electron density at mid-latitudes far more accurately 
than the simple Chapman form. We assume here that such a function 
also describes the vertical electron density distribution at the 
equator. We can then write 


f(r o ,rr/2) = N r exp V 2 {l - 


r - r 
o mo 


,r, , (r o - r mo> 

‘oL 1 ' “ exp( ' c ' 2H 



where H is the scale height of the ionic constituent and N r 
is the value of electron density at the equatorial height r mo « 
The parameter a, which is a measure of departure from the simple 
Chapman function, and thereby a shape factor, is defined as 
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a 


H - H(r ) 
o ' mo' 

— h: 


(34) 


where H(r mQ ) is that value of H(r Q ) at r Q = r mQ . Also, r Q 
is understood to be the radial height specifically at 9 = tt/2. 

Although it will not be shown here, (31) and (32) produce 
nearly identical results in the equatorial region because the 
small variation of r in the height region of our interest. 
Furthermore, the simplified form given by (31) is more convenient 
for comparison with the results of the GKS paper. We therefore 
substitute (33) into (31) and obtain 


N(r,9) = N r exp Vz -fl - 
mo 


r esc 9-r 


mo 


r o (r esc 9-r ) _ 

H^l - « exp< m )] 


o 

r cot 9 


H i 


- exp 


[- 


r esc 9-r 


mo 




r ar(r esc 9-r , _ 

2H o[ 1 ’ « ex P< 2H } ] 


(35) 


Equation (35) then provides a general expression for the electron 

density at all heights and colatitudes provided that we are 

in a region where the effects of collision can be neglected. 

The variation of N(h,9)/N. with colatitude at constant 

“mo 

height is shown in Figures 1-4 for various values of a,li ^ and 

H q . In these figures, we have converted radial height r to 

altitude h by taking the earth's radius as 6370 km. We have 

also assumed that H = H. because the effective scale height 

ox 

in (33) approaches H q at high altitudes and Chandra (1963) 
has indicated that H q becomes equal to EL at these altitudes. 
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We first note that the basic features of the theoretical 
electron density distribution are unaltered from those first 
obtained in GKS to describe the geomagnetic anomaly in the 
vicinity of the equator. Once again the theoretical description 
breaks down in the bottomside but this is precisely the region 
where the neglect of momentum transfer terms becomes in- valid. 
Furthermore, comparison of Figures 3 and 4 clearly shows the 
insensitivity of the topside results to the parameters h mQ and 
H (except for shifting the constant height profile scale 
vertically). We therefore conclude the principal properties 
of the curves can be studied quite extensively by simply altering 
the shape factor a. 

The changes due to variations in or are shown by comparison 
of Figures 1, 2 and 3. We have also provided a more detailed 
comparison for one particular height profile in Figure 5. Al- 
though we have included values up to <* = 0.6 to demonstrate 
the trend of the curves, the highest values are extreme and 
not likely to be representative of ionospheric conditions. On 
the other hand, a = 0.1 to <* = 0.4 are very reasonable values 
for us to expect under normal conditions representing diurnal 
and solar cycle variations. 

Finally, in equation (35), if we identify the term 
H o [l-a exp(- o' (r csc 2 e-r mo )/2H o ) ] with 1/k of the GKS paper, we 
see that since or is positive, kHL>l is true for all heights. 

In particular, if a = 0, we generate curves which are identical 
to those in GKS for kH ± =l. This explains why ^>1 provides 
the closest fit with experimental data in that paper. 

PROBLEMS INVOLVED IN THE DERIVATION OF THE DIFFUSION EQUATION 

In the previous sections we have seen how the equations 
of motion for electrons and ions are sufficient to obtain a 
theoretical description of the electron density distribution 
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in the topside equatorial region of the ionosphere under equinox 
conditions. This has required us to make certain assumptions 
concerning collision frequencies or velocity components along 
field lines and also forced the application of an empirical 
boundary condition at the equator. In order to produce the 
empirical boundary condition theoretically and also obtain a 
solution which is valid in both the topside and bottomside 
equatorial F region, it is necessary to turn to the continuity 
equation for additional information. Using the explicit ex- 
pressions for velocity which are derivable from the equations 
of motion, it is then possible to derive the diffusion equations 
associated with the ionosphere. 

If we simply require total ambipolar diffusion (equation 1) 
to occur in the ionosphere so that v is independent of the 
electric field explicitly, and also demand that, in all regions 
concerned, the motion along field lines are much larger than 
the drifts normal to field lines, we must then invoke v x 3 = 0 
which, using (4), (5) and (28), gives the constraint equation 


1 3N 1 2 3N 

N dr 2H N rdG 


cot 9 


(36) 


This leads to a hydrostatic distribution of electron density 
which does not agree with measured results, as has been demon- 
strated in C-I. A second approach (Kendall, 1962, and GS-II) 
is the assumption that ambipolar diffusion exists only along 
field lines (see equation (19)). Thus, if we assume that the 
parallel velocity components of electron and ion velocity are 
equal and much greater than either of the unequal perpendicular 
velocity components, we can write (20) as a good approximation 
for the entire velocity. In mathematical notation, we have 


v 





v . 
i 


X 


(37) 
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where v and v. are the perpendicular components of electron 
and ion macroscopic velocities respectively. This implies that 

V sa v (38) 

1 i 

and provides us with a velocity expression independent of electric 
field. The general contention has been that (38) allows us to 
write the steady-state continuity equation in the following form: 

Q - L = V . N v v . N v (39 

n 

where Q and L are production and loss respectively. The procedure 
has been to substitute (20) into (39) and obtain the well known 
form of the two dimensional diffusion equation without invoking 
the equation of constraint (36) . 

We wish to discuss this approach by first questioning the 
validity of (37) , and then demonstrating that even if it were 
true, (38) cannot in general imply (39) without the additional 
inclusion of the constraint equation. This will demonstrate 
that the field line ambipolar diffusion approach with neglect 
of the perpendicular velocity components is identical to the 
total ambipolar diffusion case in which velocities are assumed 
to lie along field lines. Thus, the results of the two approaches 
are identical, leading to the conclusion that ambipolar diffusion 
in which the macroscopic velocity lies along a field line, 
cannot be the correct physical model to describe the equatorial 
electron density distribution in the F region of the ionosphere. 

Let us first consider (37) . We have already seen that 
the assumption of diffusive equilibrium along a field line 
(v = o) leads to a correct description of the electron density 
in 1 at least the topside region of the ionosphere. If this is 
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the true model of the physical situation, then it is inconsistent 
with (37) and we cannot expect any results obtained using (37) 
to provide us with correct results concerning this region. If, 
on the other hand, the neglect of momentum transfer terms can 
be attributed to small collision frequencies instead of diffusive 
equilibrium, (37) need not be violated. This might be a further 
justification for validating the collision frequency assumption 
instead of the diffusive equilibrium model. Unfortunately, as 
we approach the equator, we see from (20) that v approaches 

. li 

zero since both sin I and 9N/B0 approach zero. The latter 
condition is based strictly vpon the empirical condition of symmetry 
about the equator. We therefore find that no matter how small 
v e and v^ may be, there will always be a region about the 
equator in which (37) does not apply unless 

v = v. = o (40) 

A J. 

which is identical to the equation of constraint, (36). 

We now return to the second question. That is, even if 

the parallel components of electron and ion velocities are much 

greater than the perpendicular component, which could still be 

possible provided v» ^ v* , is it possible to describe the 

li n 

electron density distributions in the entire region of the 
ionosphere by (39)? We note that 

V . Nv = v . VN + NV . v = ( v + v ) . VN + NV . (v +v ) (41) 

11 X JL 11 

where v is now either the electron or ion velocity and 

v = v + v (42) 

n J_ 

Now, in order to write (39), we must demand that 
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v . 7N + N7 . v » v . 7N + N7 . v, (43) 

11 n A. 

Although (43) could be true for certain special cases, there is 
no a priori guarantee that (43) will be implied by (37) in 
general without the additional condition that v = o. Thus, 
if we are to write (39) as a direct and general implication of 
(37) , we are once again forced to employ the constraint equation. 

We cannot state that (40) holds in a very small region 
about the equator so that its effect outside this region can 
be neglected. The geomagnetic anomaly itself is a second order 
effect and we cannot expect to reproduce it by neglecting 
second order terms which are responsible for its existence. 

We therefore find that if ambipolar diffusion exists in 
the ionosphere, and if it is restricted to the field line 
direction, we cannot assume (37) without imposing an additional 
constraint equation . Furthermore, (37) does not generally imply 
(39) in the ionosphere with or without ambipolar diffusion 
unless the constraint equation is also employed. However, since 
(37), (39), and the assumption of ambipolar diffusion along a 
field line do not provide to the correct description of equa- 
torial electron density, we must conclude that these assumptions 
are not valid in a theory leading to a description of the electron 
density distribution in the equatorial ionosphere. 

Kendall (1962), and Rishbeth, Lyon and Peart (1963), have 
attempted to numerically integrate (39) derived from (20) and 
(37), without invoking the equation of constraint. They have 
been unable to obtain the correct description of the geomagnetic 
anomaly and have therefore concluded that diffusion may not be 
a very important physical process governing the measured dis- 
tribution of electron density. However, on the basis of the 
discussion presented in this section, it now appears that the 
physical assumptions used in deriving the form of the continuity 
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equation used in their work may not be Valid, which simply 
implies that the diffusion equation is far more complicated 
than originally believed. Since (19) and (37) are no longer 
valid, we can no longer equate electron and ion velocities 
to eliminate electric field. Instead, we must write separate 
continuity equations for electrons and ions and describe the 
behavior of electric field before it is possible to obtain the 
correct theoretical description of the geomagnetic anomaly. 

It may appear that the results presented in GS-II are 
also not valid for the reasons discussed above. However, a 
closer inspection of GS-II shows that no new information was 
obtained from the solution of continuity equation than that 
already available from the equations of motion. The equation 
discussed in GS-II was simply 

V . Nv =0 (44) 

i l 

where the explicit production and loss terms were neglected in 

obtaining the series solution. Furthermore, as shown in GKS, 

the equation of motion leading to (26) , whether derived assuming 

v = 0, or by making assumptions concerning collision terms, 

has the identical solution to that obtained from (44) in GS-II. 

For the case v =0, (44) obviously cannot give any new infor- 
i i 

mat ion. This explains why the empirical boundary condition 
was necessary to obtain a non-arbitrary solution from (44) in 
GS-II. We should point out, however, that solutions of (39) 
making use of explicit production and loss terms should not 
give correct results in the equatorial regions of the ionosphere for 
the reasons discussed in this section. 

THE DISTRIBUTION OF THE NEUTRAL ATMOSPHERE 

In this section we will show that when the drag forces 
are not negligible, as might be the case in the lower F-region 
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and E-region, it is possible to study the behavior of the neutral 
atmosphere without imposing any restriction on the velocities 
of the various constituents. To obtain the necessary starting 
equation, we first sum (3), (4) and (5): 

-V(p e + p ± + P n ) + (n n m n + Nm i ) g + 3 x § = 0 (45) 

where we have once again used (14). The component of (45) along 
the direction of magnetic field is then 

[-7(p e + P ± + P n ) + (m n n n + m i N) g] . K = 0 (46) 

Comparison of (45) and (46) shows that the net force due to 
pressure gradient and gravity of all particles is perpendicular 
to the magnetic field and balanced by a current flow force. 

Next, using (7), (8), (9), (18) and (28), we have 

a (n + 2N) tan 9 9 (n + 2N) n N 

H + + _£ + — = 0 (47) 

9r r 9 9 H n EL 

where H n is the scale height of the neutral atmosphere. 

Since N « n and H. « H , we can write 
n i n' 



where the suffix on n has been dropped for simplicity. Then 
(47) becomes, in total derivative form, 


d(n .+ 2N) , n + 2N 

d? + H 

n 


Integration of (49) along the field line gives 
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n (r , 0 ) « n + 2N 


-r 

g(r 0 ,Tr/2)e r o 


dr 


(50) 


where g(r Q ,TT/2) is an arbitrary function of height at the equator 
and r Q is defined in (27) . If we now demand that the radial 
distribution of the neutrals obey the normal hydrostatic law 
at the equator, so that 


-1*0 dr 

g(r 0 ,w/2) = n co e Jr oo H n ( 51) 

where n QQ is the neutral number density at height r Qo on the 
equator, then 


_pr dr 

n(r) = n oQ e Jr oo H n, (52) 

a result which is entirely independent of 0. If, on the other 
hand, g(r Q ,n/2) is perturbed in any manner from the exact hydro- 
static equilibrium case, we will obtain a distribution for n 
which does depend on 0. The origin of this angular dependence 
on the neutrals may seem somewhat surprising until we realize 
that in selecting a functional form for g(r o ,Tr/2), any deviation 
in the equatorial neutral distribution from hydrostatic equi- 
librium must arise due to collisions between neutrals and geo- 
magnet ically controlled charged particles. Thus, if the col- 
lisions between neutrals and charged particles are sufficiently 
large to make the momentum transfer forces between charged and 
neutral particles important, the neutrals will begin to tend 
toward the angular distribution of the geomagnet ically controlled 
particles. This can also be seen from (3), where it is obvious 
that we will not obtain the exact hydrostatic distribution in 
a region when the terms on the left hand side become important. 


20 


On this basis, we might expect to observe angular variations 
of the neutral distribution in the bottomside regions of the 
ionosphere where charged-neutral particle interactions become 
important . 

CONCLUSIONS 

From the discussion and results of this paper, we have 
shown the following: 

1. From the equations of motion, it is possible to derive 
an expression for the electron-density distribution along a 
field line either by assuming diffusive equilibrium along the 
direction of the magnetic field or by neglecting the drag forces 
arising from collisions. The latter assumption appears to be 
more realistic in the topside of the ionosphere. In either case, 
it is necessary to assume a radial distribution at the equator 

to obtain the electron density distribution. 

2. We have provided a more accurate formula for the repre- 
sentation of the equinox geomagnetic anomaly than that produced 
in GKS. Since the empirical boundary condition equation used 
herein has been shown by Chandra (1963) to fit nearly all vertical 
profiles of electron density measured to date, we can safely 
assume that the proper selection of parameters in this formula 
will lead to a reasonable reproduction of the anomaly in any 
equatorial region of the ionosphere where interactions of neutrals 
with charged particles are small because of infrequent collisions. 

3. The theory discussed above is semiphenomenological ; 
i.e., it is based on effect and not cause. It does not require 
a knowledge of the complicated array of physical effects and 
mechanisms which combine to form the geomagnetic anomaly but, 
instead, uses an empirical boundary condition which is the 
accumulated effect of all these causes. 

Naturally, if we are to increase our knowledge of 
the basic mechanisms causing the anomaly and thereby replace 



21 


the empirical boundary condition by one based on more funda- 
mental considerations than measurement, we must turn to the 
equations of continuity. Unfortunately, the derivation of 
the correct continuity equations requires knowledge concerning 
the electron and ion velocities and/or the electric fields acting 
on these particles. Currently, most derivations of the conti- 
nuity equations employ simplifying assumptions, such as 

v = v. = v n ! v n » v e > v - • The equation derived in 
®xi Hj e X *x 

the literature under the above assumptions has been numerically 

integrated by several workers to obtain a theoretical electron 
density distribution near the equator under steady state condi- 
tions. The results obtained by these workers have been unable 
to account for the gross features of the geomagnetic anomaly, 
at least to the correct order of magnitude. This has lead them 
to believe that diffusion is of minor importance in governing 
the geomagnetic anomaly. 

We have been able to demonstrate that the velocity assump- 
tions described above do not lead to the proper description 
of the geomagnetic anomaly. We therefore feel that the assump- 
tion about velocities used in the continuity equation rather 
than the ineffectiveness of motions are responsible for the 
unsatisfactory description of the geomagnetic anomaly obtained 
by others. 

4. A study of the neutral atmosphere distribution has 
led us to the conclusion that geomagnetic control of neutrals 
occurs in any region of the ionosphere where interactions of 
neutrals with charged particles become important. Since this 
is most likely to occur in the lower F region of the ionosphere 
we suggest that such geomagnetic control of the neutrals might 
be observable in this region. 
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LIST OF FIGURES 

Figure 1 - Constant Height Profiles of Relative Electron Density 

vs. Colatitude for a = 0, h,_ = 500 km, H = 100 km. 
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Figure 2 - Constant Height Profiles of Relative Electron Density 

vs. Colatitude for oi = 0.1, h_ - 500 km, H = 100 km. 
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Figure 3 - Constant Height Profiles of Relative Electron Density 

vs. Colatitude for a = 0.3, h = 500 km, H = 100 km. 
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Figure 4 - Constant Height Profiles of Relative Electron Density 

vs. Colatitude for a » 0.3, li „ = 400 km, H = 75 km. 
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Figure 5 - The cn Dependent Behavior of the 500 km Profile of 

Relative Electron Density vs. Colatitude for 

h _ a* 500 km, H = 75 km. 
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